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Landau hydrodynamics is a plausible description for the evolution of the dense hot matter pro- 
duced in high-energy heavy-ion collisions. We review the formulation of Landau hydrodynamics 
to pave the way for its application in high-energy heavy-ion collisions. It is found that Landau's 
rapidity distribution needs to be modified to provide a better quantitative description. In par- 
ticular, the rapidity distribution in the center-of-mass system should be more appropriately given 
^q" ■ as dN/dy oc exp{^/y^ — y 2 }, where j/j = ln{^/sjvjv /m p } is the beam nucleon rapidity, instead of 

Landau's original result of dN/dy (Landau) oc exp{ \J L 2 — y 2 } where L = lnl^/sjvjv /2m p }. The 
modified distribution is compared with the Landau distribution and experimental data. It is found 
£N| ' that the modified distribution agrees better with experimental dN/dy data than the Landau distribu- 

tion and it differs only slightly from the Landau Gaussian distribution dN / 'dy(Landau Gaussian) oc 
exp{— y 2 /2L}. Past successes of the Gaussian distribution in explaining experimental rapidity data 
^yy • arises, not because it is an approximation of the original Landau distribution, but because it is in 

fact a close representation of the modified distribution. 
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I. INTRODUCTION 

6 
=5 

The evolution of the dense hot matter produced in high-energy heavy-ion collisions is a problem of 
intrinsic interest. In many problems in high-energy heavy-ion collisions, such as in the investigation of 
the fate of a heavy quarkonium in quark-gluon plasma and the interaction of a jet with the produced 
t*** . medium, it is desirable to have a good description of the dynamics of the matter produced in the 
collision. For these problems, Landau hydrodynamics [l|, 0] furnishes a plausible description for the 
evolution of an assembly of dense matter at a high temperature and pressure. Its dynamics during the 
ON ■ first stage of one-dimensional longitudinal expansion can be solved exactly @, 1,1, 0, 0, 0, 1, §, 0, 
The one-dimensional longitudinal expansion problem admits an approximate solution that 



q \ is applicable to the bulk part of the fluid. The subsequent three-dimensional motion can be solved 
approximately to give rise to predictions that come close to experimental data 13, 0, 15, 16] . A 



thorough understanding and critical tests of the model will make it a useful tool for the description 
of the evolution of the produced dense matter. We therefore review the formulation of Landau 
c$ \ hydrodynamics to pave the way for its quantitative applications in high-energy heavy-ion collisions. 
Common to both Landau hydrodynamics and Bjorken hydrodynamics is the basic assumption 
that in the dense hot matter produced in high-energy heavy-ion collisions, the density of the quanta 
of the medium (either quarks and gluons or the hadrons) is so high that the mean-free paths of 
the constituents are short, and a state of local thermal equilibrium can be maintained through out. 
Such a condition may be at hand in collisions at very high energies. Under the circumstances, 
the dynamics of the system will be determined by the equation of hydrodynamics possessing local 
thermal equilibrium. As we are dealing with high temperatures and high densities for which the 
temperature is of the order of the rest mass of the quanta, we need to examine the problem with 
relativistic hydrodynamics. 
The major difference between Landau hydrodynamics and Bjorken hydrodynamics is in the nature 
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of the flow. In Bjorken hydrodynamics, boost invariance is assumed so that the energy density e 
is independent of the rapidity. It is a function of the proper time r only, eBjorkcnfV)- However, in 
Landau hydrodynamics, no such an assumption is made, and the energy density is a local function 
of spatial and temporal coordinates, eL an dau(£, z, x, y). 

We can consider collisions of two equal-size nuclei of diameter a. The disk of initial configuration 
in the center-of-mass system has a longitudinal thickness A given by 



A = a/ 7 , 



(i-i; 



where 7 is the Lorentz contraction factor 
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2m r 



;i.2) 



y/s NN /2 is the center-of-mass energy per nucleon, and m p is the proton mass. The dimensions of 
the disk are not symmetrical in the longitudinal and transverse directions. The initial configuration 
can be depicted as in Fig. 1 where it can be represented as an elliptical disk with a thickness A and 
major diameters a x and a y , with the reaction plane lying on the x-z plane. Depending on the impact 
parameter, the dimensions of the disk obey a x < a y < a, and the azimuthal radius depends on 
the azimuthal angle <p measured from the x-axis. For a central collision, a x = a y = a. 




Initial configuration 



FIG. 1: Initial configuration in the collision of two heavy nuclei in the center-of-mass system. The region of nuclear overlap 
consists of a thin disk of thickness A along the longitudinal z-axis. The reaction plane is designated to lie on the x-z plane, 
and the transverse radii are a x /2 and a y /2. 



The matter in the disk can be described as a fluid with an energy density field e, and a 4-velocity 
field u^. In the center-of-mass system to which we focus our attention, the disk of matter can be 
considered to be initially nearly at rest. The density of the produced matter and its associated 
pressure are very large, whereas the spatial dimensions of the fluid are very small. As a consequence, 
large pressure gradients are generated. The generated pressure gradients lead to the subsequent 
expansion of the system in all directions. 

In highly relativistic collisions to which we focus our attention, the center-of-mass energy per 
nucleon is much greater than the nucleon rest mass, (a/stvat/2) m p . We have then A -C a x < a y 
and V z p ^> V-rP > V^p. The pressure gradient in the longitudinal z direction far exceeds those 
along the transverse x and y directions. The expansion along the longitudinal direction proceeds 
much faster than the transverse expansions. Among the two transverse directions, the expansion 
along the x direction proceeds slightly faster than the expansion along the y direction, leading to 
the presence of the elliptic flow. 
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The large difference of the expansion speeds in the longitudinal and transverse directions al- 
lows Landau to split the evolution dynamics into two stages. In the first stage, one considers a 
one-dimensional longitudinal expansion in conjunction with a simultaneous but slower transverse 
expansion. The second stage occurs when the transverse displacement exceeds the initial trans- 
verse dimensions so that the transverse force acting on the fluid element become so small as to be 
negligible. Thus, in the second stage, the fluid element will freeze its rapidity to stream out in a 
conic flight with a fixed polar angle. The magnitude of the transverse displacement decreases with 
increasing rapidity magnitudes. Therefore different parts of the fluid will enter from the first stage 
to the second stage at different times. Those fluid elements with the smallest rapidity magnitude 
will enter the second stage the earliest while those elements with the largest rapidity magnitude will 
enter the second stage the latest. By matching the solutions for the first and the second stages, the 
dynamics of the system can be followed, up to the end point of the hydrodynamical evolution. 

Recent comparison with experimental data indicates that the Landau hydrodynamical model yields 
results that agree with experiment. In Refs. Il3|, [14|, [l5|, [l7| , quantitative analyses use an approximate 
form of the Landau rapidity distribution [Hill], 



dN/dy oc exp{-y 2 /2L}, 
where L is the logarithm of the Lorentz contraction factor 7 

L = In 7 = \n(y/s NN /2m p ) 



;i.3) 



;i-4) 



This Gaussian form of the rapidity distribution gives theoretical rapidity widths that agree with 
experimental widths for many different particles in central AuAu collisions, to within 5 to 10%, 



from AGS energies to RHIC energies [13|, [14|, [15j . The Landau hydrodynamical model also gives the 



correct energy dependence of the observed total multiplicity, and it exhibits the property of "limiting 



fragmentation" at forward rapidities [3, lEf. A similar analysis in terms of the pseudorapidity 
variable 77 at zero pseudorapidity has been carried out in [l8| . 

The successes of these quantitative and qualitative analyses indicate that Landau hydrodynamics 
can be a reasonable description. However, they also raise many unanswered questions. Firstly, the 
original Landau result stipulates the rapidity distribution to be [H, 0] 



dN/dX (Landau) oc expjV L? — A 2 }, 



;i.5) 



where the symbol A is often taken to be the rapidity variable y [HI, 14, 15, T3]. In the original 
work of Landau and his collaborator in [H, 0] , the variable A is used to represent the polar angle 9 as 
e~ A = 9; there is the question whether the variable A in the Landau rapidity distribution (jl.5p should 



be taken as the rapidity variable y [13|, [14|, [15|, [17| or the pseudorapdity variable 77 [18] appropriate 



to describe the polar angle. Such a distinction between the rapidity and pseudorapidity variables is 
quantitatively important because the shape of the distributions in these two variables are different 
near the region of small rapidities 19| . Secondly, the Gaussian rapidity distribution (II. 3p used in the 



analyses of Refs. [l3j, ll4|, [l5j], as well as in the earlier work of (l7j|, is only an approximation of the 
original Landau distribution of Eq. (11.51) . The Gaussian approximation is valid only for A C L, when 
we expand the square- root function in powers of X/L. There is the question whether it is appropriate 
to apply the Gaussian distribution (I1.3P to the whole range of A, including the region of A for which 
the expansion condition A <C L is not satisfied. The Landau distribution (II. 5p and the Gaussian 
distribution (II. 3p are in fact different distributions. While the original Landau distribution may be 
considered to receive theoretical support in Landau hydrodynamics as justified in Refs. [H, 0], a firm 
theoretical foundation for the Gaussian distribution (jl.3p in Landau hydrodynamics is still lacking. 
Finally, if one does not make the expansion of the square-root function and one keeps the original 
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Landau distribution of Eq. (]1.5p . then there is the quantitative question 2(| whether this original 
Landau distribution with the theoretical L value will give results that agree with experimental data. 

In view of the above unanswered questions, our task in reviewing the Landau hydrodynamical 
model will need to ensure that we are dealing with the rapidity variable y and not the pseudorapidity 
variable. We need to be careful about various numerical factors so as to obtain a quantitative 
determination of the parameters in the final theoretical results. Finally, we need to ascertain whether 
the theoretical results agree with experimental data. If we succeed in resolving these questions, we 
will pave the way for its future application to other problems in high-energy heavy-ion collisions. 



II. TOTAL NUMBER OF PRODUCED PARTICLES 



Before we study the dynamics in detail, we would like to examine the total number of particles 
produced in a central heavy-ion collision in Landau hydrodynamics. Landau assumed that the hy- 
drodynamical motion of the fluid after the initial collision process is adiabatic. He argued that the 
only thing that can destroy adiabaticity would be the shock wave. They occur at the initial com- 
pressional stage of the collision process [2lj]. It is hard to imagine how they could be formed in the 
bulk part during the subsequent expansion phase after the initial compression and thermalization. 
Landau therefore assumed that during the longitudinal and transverse expansion phase under con- 
sideration, the entropy content of the the individual region remains unchanged. The total entropy 
of the system is therefore unchanged and can be evaluated at the initial stage of the overlapped and 
compressed system. 

The total entropy content is a very useful quantity because it is closely related to the total particle 
number. From the consideration of the thermodynamical properties of many elementary systems, 
Landau found that the ratio of the entropy density to the number density for a thermally equilibrated 
system is nearly a constant within the temperature region of interest. Landau therefore postulated 
that the number density is proportional to the entropy density. Thus, by collecting all fluid elements, 
the total number of particles is proportional to the total entropy. As the total entropy of the system 
is unchanged during the hydrodynamical evolution, the total number of observed particles can be 
determined from the initial entropy of the system. 

We work in the center-of-mass system and consider the central collision of two equal nuclei, each of 
mass number A, at a nucleon-nucleon center-of-mass energy y/s NN . Consider first the case of central 
AA collisions with A ^> 1 such that nucleons of one nucleus collide with a large numbers of nucleons 
of the other nucleus and the whole energy content is used in particle production. This is the case of 
"full stopping" . The total energy content of the system is 

E = ^A. (2.1) 

The initial compressed system is contained in a volume that is Lorentz contracted to 

Air 

V = -{a/2fh, (2.2) 

where the nuclear radius a/2 is related to the mass number by 

a/2 = roA 1 / 3 , (2.3) 
and r = 1.2 fm. The energy density of the system is therefore 



6 = E/V = 7v ^7/(47rr 3 /3), 



(2.4) 



which is independent of A and depends on energy as s NN . For a system in local thermal equilibrium, 
the entropy density a is related to the energy density by 



cr = constant e 3//4 . 
The total entropy content of the system is therefore 

S = aV = constant s 1//4 A. 



(2.5) 



(2.6) 



With Landau's assumption relating entropy and particle number, N oc S, the total number of 
particles produced is 



N oc s^A, 

and the total number of produced charged particles per participant pair is 

N ch /A = iV ch /(iV part /2) = Ki^/GeV) 1 / 2 , 
where K can be determined phenomenologically by comparison with experimental data. 

40 



(2.7) 
(2.8) 
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FIG. 2: Total number of produced charged particles per pair of participants, A r c h/(A ? par t/2), as a function of (^/s JVJV /GeV) 1//2 . 
(a) PHOBOS data A r c h/(A r par t/2) data for central AuAu collisions at different (y/s NN /GeV) 1 ^ 2 and the Landau hydrodynamical 
model fit, and (b) the extrapolation of the Landau hydrodynamical model to LHC energies. 

In Fig. 1(a), we show the PHOBOS data of A r c h/(A r part/2) as a function of (v^^/GeV) 1 / 2 for 
central AuAu collisions in RHIC 

13 EE IS- 

The RHIC AuAu data can be parametrized as 



AW(^Wt/2) = 1.135 + 2.019( v ^7/GeV) 1 / 2 , 



(2.9) 



where the constant 1.135 arise from the leading nucleons. The constant K as determined from the 
data is K = 2.018 which agrees with the earlier estimate of K = 2 [l], 0]. 
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Consider next pp and pp collisions in which not all the energy of y/s NN is used in particle pro- 
duction, as the leading particles carry a substantial fraction of the initial energy. If we denote the 
particle production energy fraction in pp and pp collisions by £, then Eq. (I2.8P is modified to be 

= ^(eV^/GeV) 1 / 2 . (2.10) 

Comparison of the charged particle multiplicity in pp, pp, and e + e~ collisions indicates that the 
particle production energy fraction £ for pp and pp collisions is approximately 0.5 [l4, 15, 19, 0. In 



contrast, the case of RHIC AA data in high-energy heavy-ion collisions corresponds to full nuclear 
stopping with £ = 1 [3, IEI . 



In Fig. 1(b), we show the predictions for the charge particle multiplicity per pair of participants 
for collisions at LHC energies. For pp collisions at 14 TeV with a particle production energy fraction 
£ = 0.5, iV c h is predicted to be 170. For central PbPb collisions at y/s NN = 5.5 TeV with full nuclear 
stopping (£ = 1), A r ch /(A r part /2) is predicted to be 151. 



III. LONGITUDINAL HYDRO DYNAMICAL EXPANSION 



We proceed to examine the dynamics of the longitudinal and transverse expansions. Among the 
coordinates (t, z, x, y) = ( X - X - X - X ) used as the arena to describe the fluid, Landau suggested a 
method to split the problem into two stages. The first stage consists of simultaneous and independent 
expansions along the longitudinal and the transverse directions. For the one-dimensional longitudinal 
expansion, the equation of hydrodynamics is 

dT oo dT oi 

where 

= (e + p)uV -pg^. (3.3) 

In hydrodynamical calculations, we need an equation of state. We shall assume for simplicity the 
relativistic equation of state 

P = e/3, (3.4) 

which allows one to express the local pressure p as a function of the energy density e. It is convenient 
to represent the velocity fields (u°, u 1 ) by the flow rapidity y 

u° = cosh?/, (3.5a) 
u = siiahy. (3.5b) 

In this form, the property of the velocity fields, (u ) 2 — (u 1 ) 2 = 1, is readily satisfied. A par- 
ticle of mass m flowing with such a rapidity y will have a momentum (p ,^ 1 ) = (mvP,mu l ) = 
(m cosh y, m sinh y) . 
It is convenient to introduce the light-cone coordinates t + and t- 

t+ = t + z, (3.6a) 
t_ = t-z, (3.6b) 
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and their logarithmic representations (y + , yJ) defined by 

y ± = \n{t ± /A} = \n{(t±z)/A}. (3.7) 

The advantage of using the these variables is that the hydrodynamical equations and the corre- 
sponding solutions are quite simple in these variables. The knowledge of the fluid energy density e 
and the flow rapidity y, as a function of (t + , t_) or (y + , y_) will provide a complete hydrodynamical 
description of the whole system. 
Exercise (1.1) shows that the hydrodynamical equations (13.11) and (13. 2p become 

de d(ee~ 2y ) 



dt+ dt 
d(ee 2y ) de 



dt, dt. 



0. (3.8b) 



The t + and t- coordinates are light-like because they lie on the light cone. If we view the in- 
coordinate as approaching time-like and the t_ as approaching space-like, then the first equation can 
be considered the equation of continuity and the second the Euler equation of the hydrodynamical 
current. Conversely, if we view the t_ coordinate as approaching time-like and the t + as approaching 
space-like, then the second equation can be considered the equation of continuity and the first the 
Euler equation of the hydrodynamical current. 

For the first stage of one-dimensional hydrodynamics, the exact solution for an initially uniform 
slab has been obtained [s| and discussed in 0, y, 0, [H|, @, Q, H, M, EU U, 12]. There is in addition an 
approximate particular solution that is very simple and applicable to the bulk part of the fluid [H, 0]. 
In view of the matching of the solution to an approximate three-dimensional motion in the second 
stage, it suffices to consider the approximate solution of longitudinal hydrodynamical expansion in 
the present discussions. 

The approximate particular solution of the hydrodynamical equation is 0] 

t(y+,v-) = e exp + y-- Vv+y^)^ , (3.9a) 

y(y + ,y-) = (y+-y_)/2. (3.9b) 
The rapidity solution can also be written alternatively as 

e 2y(y + ,y-) = U = t±l_ ( 3 . 10 ) 
t_ t-z y ' 

The corresponding velocity field u° is related to the fluid coordinates by 

u\y + ^) = -=S 7 == . (3.11) 

y/t - Zy/t + Z 

The constant eo in Eq. (I3.9al) is related to the initial energy density at (y+o, y~o) by 

eo = e(2/ + o,y-o)e^, (3.12) 

where O is 

4 

4>o = g(y+o + y-o - Vy+oy-o)- ( 3 - 13 ) 
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We can easily prove by direct substitution that (I3.9ap and f!3.9bj) (or (I3.10p ) are approximate par- 
ticular solutions of the hydrodynamical equations (I3.8al) and ( I3.8bl) . First, substituting Eq. (13.101) 
into the hydrodynamical equations, we obtain 



9f „ 

at + +I 


" de 

_dt. 


i " 


t_ 
U 




" de 


i " 


4 + + 


de 


2 


_dt + + 


*+. 




at. 



0. 
0. 



(3.14a) 
(3.14b) 



We write out t_/t + in the second equation and substitute it into the first equation, and we get 



de de 
dt + dt- 



4 



de_ 
dL 



de_ 

dtl 



1 







(3.15) 



We multiply this expression by and change into the logarithm variables y + and ?/_, then the 
above equation becomes 



de de 
dy+ dy^ 



1 



de 
dy~ 



1 



de 
dy+ 







(3.16) 



If we now substitute Eq. (I3.9ap for e into the lefthand side of the above equation, we find that the 
lefthand side gives zero, indicating that Eqs. (I3.9al) and (I3.9b|) are indeed the approximate particular 
solution of the hydrodynamical equation. 

The above solutions (I3.9ap and (13.9bp are applicable for both positive and negative z values and 
contains the proper reflectional symmetry with respect to the interchange of y + and The so- 
lutions are appropriate for the description of the dynamics of the dense hot matter in the first 
stage of longitudinal expansion. For simplicity, we shall concentrate our attention on one of the two 
symmetric sides, the side with positive z. 

The simple approximate solution of f)3.9ap and (13.9bp have limitations. Because of the nature of 
the approximate solution, it cannot describe the boundary layer for which t — z < A and y_ becomes 
negative. At the boundary layer, thermodynamical quantities such as the energy density and the 
entropy density decrease rapidly and it is not included into the region for the applications of the 
solution of Eqs. (I3.9al) . In highly relativistic collisions A/a is given by (2m p / ■ s /s~^), which is a 
small quantity. The region excluded from the approximate solution is not significant in a general 
description of the fluid. 

There is another minor disadvantage of the solution. Because of the nature of the solution for 
the bulk matter, the solution in Eq. f)3.9ap provides only limited choice on the initial conditions, 
within the form as specified by the simple functions in these equations. However, the dynamics of 
the system is likely to be insensitive to the fine characteristics of the initial density distribution of 
the contracted disk of thickness A along the longitudinal direction. A thin slab of matter of the 
right dimensions within the Landau model will likely capture the dominant features of the evolution 
dynamics. 

The solutions e(t + , t_) and y(t + , t_) can be transformed into e(t, z) and y(t, z) by a simple change 
of variables. We show the solution e(t, z) and y(t, z) as a function of z/A for different t in Fig. 2(a) 
and 2(b) respectively. As one observes, the maximum energy density decreases nearly inversely as a 
function of t. One notes that the energy density in the central region decreases in time approximately 
as 1/t, as the matter expands outward. The flow rapidity is initially close to zero. As the bulk matter 
expands longitudinally, the magnitude of the rapidity field increases at large values of |z|/A. We 
show the velocity field u°(t,z) and the corresponding longitudinal velocity v z (t,z) as a function of 
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FIG. 3: Approximate Landau hydrodynamical solution of the one-dimensional longitudinal expansion of the system with initial 
energy density to, as a function of the longitudinal coordinate z and various times t in unit of the slab thickness A. 
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FIG. 4: The proper velocity u and the velocity v z along the longitudinal z-direction as a function of z and t Landau 
hydrodynamics. 



z/A for different t in Fig. 3(a) and 3(b). The magnitude of the longitudinal velocity of the fluid is 
large at the region of \z\ close to t. 

The hydrodynamical solution Eq. f)3.9ap can be equivalently expressed as a function of any other 
two independent variables. In our problem, it is desirable to use the variables t and y because 
the condition for the transverse displacement to exceed the initial transverse dimension is most 
conveniently written for a fixed time t, and the dynamics after the transversality condition is satisfied 
leads to the freezing of the rapidity y. 
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In terms of the (t, y) variables, we can express y± as 

y± = ln(t/coshyA) ±y. (3.17) 
Therefore the hydrodynamical solution Eq. (I3.9ap can be re-written in terms of (t, y) as 

e(t,y) = e exp j-^(2 In (t/cosh y A) - y/[]n (t/coshy A)] 2 - y 2 ) j , (3.18a) 

z(t,y) = t sinhy/coshy. (3.18b) 

We can further express the solution as a function of (r, y) to compare Landau hydrodynamics with 
Bjorken hydrodynamics. We make the transformation 

t = r cosh y, (3.19a) 

z = rsinhy. (3.19b) 

We then have 

y± = ln(r/A)±y. (3.20) 
The energy density in Landau hydrodynamics, expressed in terms of r and y is then 

e(r,y) = e exp |~ [21n(r/A) - vX^/A)] 2 - y 2 ] | . (3.21) 

In the region y <C ln(r/A), we have 

e(r,y) = e exp |-|ln(r/A)J oc (3.22) 

which is the Bjorken hydrodynamics results. Therefore, in the region of small rapidities with \y\ <^ 
ln(r/A), Landau hydrodynamics and Bjorken hydrodynamics coincide. However, in other regions 
when \y\ is not very smaller in comparison with ln(r/A), the two types of hydrodynamics differ. In 
general, because Landau hydrodynamics covers a wider range of rapidities which may not be small, 
it is a more realistic description for the evolution of the hydrodynamical system. 

(|)[ Exercise 1.1 Show that the hydrodynamical equations for e and y are given by Eqs. (pT8a)l and IpPfoJl , 
The energy momentum tensor T 00 is equal to 

rriOO I , \ 

1 = (e + p)u u — p 

— (e + p)cosh y — p. (3.23) 

Upon using the equation of state p = e/3 in Eq. (|3.4j> . the energy momentum tensor can be expressed in terms of e and y. 

T 00 = e(cosh 2 y + i sinhy). (3.24) 

We have also, 

I = (e + PJti w 

= (4/3)ecoshysinh?/, (3.25) 

T 11 = (e + p)^ 1 ^ 1 + p 

= e(sinh 2 y + i cosh 2 y). (3.26) 
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The transformation of (t,z) to (t+,t-) lead to the following relation between the derivatives, 



a 


dt+ 


d 


dt. 


d 


d 


d 

+ dE 


m 


dt 


dt+ 


+ 1H 


dt- 


~ dt+ 


d 

dz 


_ dt+ 

dz 


d 
dt+ 


dt. 
dz 


d 
dt. 


d 

~ dt+ 


d 

dt. 



Therefore, Eqs. (pITj) and become 
dT oo dT oi g 



dt dz dt + 

dT oi dT u Q 
+ 



sinh y cosh y + e(cosh 2 y + i sinh 2 y) 



4 2 1 2 

— e sinh y cosh y + e (sinh y + — cosh y) 



+ 



at. 



dt dz dt 
By forming the difference of the above two equations, we get 



+ 



dt- 



- sinh y cosh y + e(cosh 2 y + — sinh 2 y) 



4 2 1 2 

— e sinh y cosh y — e (sinh y + — cosh y) 



dt 

This can be simplified to be 



d 1 d 8 4 

— — [e(cosh 2 y+ - sinh 2 y)] + — — [--e sinh y cosh y + -e(cosh 2 y + sinh 2 y)] = 0. 



dt. 1 3 



du +z dt- - u ' 



which is Eq. (|3.8a[) . Similarly, by forming the sum of the two equations, we get 

d 



d 
dt+ 



e(cosh 2 y + sinh 2 ) + — e sinh y cosh y 



+ 



dt- 



-e(cosh 2 y — sinh 2 y 



This leads to 



which is Eq. ([3~8b)l . 



r, d(ee 2 *) de 
dt+ dt- 



(3.27) 
(3.28) 

(3.29) 
(3.30) 

(3.31) 

(3.32) 

(3.33) 

(3.34) 
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IV. TRANSVERSE EXPANSION 



The initial configuration is much thinner in the longitudinal direction than in the transverse 
directions. Therefore, in the first stage of evolution during the fast one-dimensional longitudinal 
expansion, there is a simultaneous but slower transverse expansion. The difference in the expansion 
speeds allows Landau to treat the longitudinal and transverse dynamics as independent expansions. 
The rate of transverse expansion can then be obtained to provide an approximate description of the 
dynamics of the system. 

We shall consider first the case of a central collision, for which The case of non-central 

collisions will be discussed in Section IX. The transverse expansion is governed by the Euler equation 
along one of the transverse directions, which can be taken to be along the x direction, 



dT 02 dT 22 

+ ^— = 0, 



Of 



dx 



where 



T 02 = (e + p)u°u 2 = ^e(u°) V, 



(4.1) 



(4.2) 
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and we have used the relation u 2 = u°v x . The energy- momentum tensor T 22 is 

T 22 = (e + p){u 2 ) 2 -pg 22 = \e{u«) 2 {v x ) 2 + P- (4-3) 

As the transverse expansion is relatively slow, we can neglect the first term on the righthand side of 
the above expression and keep only the pressure term p. 

In Landau's method of splitting the equations, one makes the approximation that during the 
first stage the quantities e and y as a function of t and z have been independently determined in 
the one-dimensional longitudinal motion. It suffices to use Eq. ( 14. ip to determine the transverse 
displacement as a function of time. Equation (14. ip can therefore be approximated as 

The transverse displacement x(t) (relative to zero displacement) as a function of time t is related to 
the acceleration dv z /dt by 

The pressure is p = e/3 at the center of the transverse region and is zero at the radial surface a/2. 
Therefore the equation for the displacement is given from Eq. (14. 4p by 

4 o 2 2x^t) = e 

3 v ' t 2 3a/2 v ' 

We note that there is a factor of 4 arising from the ratio of 4e/3 from (e + p) on the lefthand side 
relative to e/3 from the pressure p on the righthand side. However, in the original formulation 
of Landau [H, 0], this factor of 4 is taken to be unity for an order of magnitude estimate of the 
transverse displacement. For our purpose of making quantitative comparison with experimental 
data, this factor of 4 cannot be neglected. We shall find in Section 5 and 6 that the presence of this 
factor of 4 leads to the modified rapidity distribution (I6.12p that agrees better with experimental 
dN/dy data than the Landau distribution of Eq. (11.51) . 

From Eq. (14. 6p . the transverse displacement x(t) during the one-dimensional longitudinal expansion 
increases quadratically with t as 

t 2 t 2 

x(t) = - . — = a-. (4.7) 

w 4a(w ) 2 4a cosh 2 ?/ V ; 

For a fixed t, the transverse displacement x(t) as a function of the rapidity is greatest at y = and 
it decreases as \y\ increases. The variation of the transverse displacement as a function of rapidity 
leads to different times when the fluid element enters from the first stage to the second stage of 
frozen rapidities. 



V. SECOND STAGE OF CONIC FLIGHT 



The description of one-dimensional solution of longitudinal expansion with the accompanying 
transverse expansion in the first stage is applicable as long as the transverse displacement x(t) is 
sufficiently small compared with the transverse dimension a. It ceases to be applicable at the time 
t m when the transverse displacement x{t m ) is equal to the transverse dimension a. Landau suggested 
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that when x(t m ) is equal to a we need to switch to a new type of solution in the second stage of 
fluid dynamics. 

What type of fluid motion is expected after the fluid element enters the second stage? For a 
fluid element at t > t m , the transverse displacement x(t) already has exceeded the initial transverse 
dimension. With the fluid element beyond the initial transverse dimension, the derivatives of hydro- 
dynamical quantities with respect to both the transverse coordinates and t_ are small. Because of 
the smallness of these derivatives, Landau argued that the hydrodynamical forces become so small 
that they can be neglected in the hydrodynamical equations at these locations. The flow rapidity y 
will be frozen for t >t m . Freezing the rapidity of a fluid element as a function of time is equivalent 
to freezing the opening polar angle 9 between the fluid trajectory and the collision axis at these 
longitudinal locations. The motion of the fluid element with a fixed polar angle can be described as 
a 'three-dimensional' conic flight, with no change of the flow rapidity. 

In mathematical terms, Landau's condition for rapidity freezeout occurs at the time t m (y) when 
the transverse displacement satisfies [H, 0] 

x(t m ) = a, (5.1) 
which, as determined from Eqs. (14.71) and (15. ip . takes place at 

tm{y) — 2au° = 2a cosh y. (5.2) 
The switching time t m (y) increases with increasing magnitude of the rapidity. 




FIG. 5: Schematic picture of the opening polar angle 9 between the fluid element trajectory and the longitudinal axis. The 
conic energy and entropy within the angle element d6 is preserved in this second stage of conic flight, 

The solution of the hydrodynamical variables during this second stage of conic flight can be easily 
obtained. In a conic flight with a cone of opening polar angle 9 within an angle element d9 as 
shown in Fig. 5, the energy- momentum tensor and the entropy flux within the cone element must 
be conserved as a function of time. The cross sectional area of such a cone element is 2irxdx. So the 
conservation of energy and entropy conic flow correspond to 

dE = eu°u°27rxdx = constant, (5.3) 

and 

dS = cru°27rxdx = e 3 ^ 4 u°27rxdx = constant, (5.4) 
where a is the entropy density. Dividing the first equation by the second equation, we get 

e 1 / 4 ^ = constant, (5.5) 

which gives 

e a (5 - 6) 

On the other hand, in the conic flight, x and dx are proportional to t. Hence, Eq. ( 15.31) gives 

eu°u°t 2 = constant. (5.7) 
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Eqs. ( 15. 6p and (15. 7p yield the dependence of various quantities as a function of t, 

e oc — , a oc — , and u° oc t. (5.8) 

These equations give the solution of the evolution of the fluid elements as a function of time in the 
second stage. By matching the solutions at t — t m (y), the energy density and velocity fields at the 
second stage for t > t m (y) is 

e(t,y) = e(t m} y)t 4 Jt 4 (5.9a) 
u°(t,y) = u°(t m ,y)t/t m . (5.9b) 

The velocity component u° of a fluid element increases linearly with time in this second stage. In 
contrast, in the first stage, the velocity component u° of a fluid element with fixed t — z and z 
comparable to t increase with time only approximately as t 1 ^ 2 , as indicated by Eq. (13.111) . Thus, the 
fluid element approaches the speed of light faster in the second stage than in the first stage. 



VI. RAPIDITY DISTRIBUTION IN HIGH ENERGY HEAVY-ION COLLISIONS 



The picture that emerges from Landau hydrodynamics can be summarized as follows. For an initial 
configuration of a thin disk of dense matter at a high temperature and pressure, the first stage of 
the motion is a one-dimensional longitudinal expansion with a simultaneous transverse expansion. 
The transverse expansion lead to a transverse displacement. When the magnitude of the transverse 
displacement exceeds the initial transverse dimension, forces acting on the fluid element becomes 
small and the fluid elements will proceed to the second stage of conic flight with a frozen rapidity. 
As the transverse displacement depends on rapidity, and the transverse displacement magnitude 
increases with increasing rapidity magnitude, the moment when the fluid element switches from the 
first stage to the second stage depends on the rapidity. The final rapidity distribution of particles is 
therefore given by the rapidity distribution of the particles at the matching time t m (y). 

We shall first evaluate the entropy distribution as a function of rapidity y and time t in the first 
stage of hydrodynamics. Consider a slab element dz at z at a fixed time t. The entropy within the 
slab element is 

dS = au°dz. (6.1) 

Using the solution (13.101) . we can express z as a function of t and rapidity y during the one- 
dimensional longitudinal expansion, 

z = t sinhy/coshy. (6-2) 

For a fixed value of t, we therefore obtain 

dS = au°tdy / cosh 2 y 

= atdy /coshy. (6.3) 

The entropy density a is related to e by a = ce 3//4 and e is given by (I3.9al) . We obtain the entropy 
element at the time t, 



dS = ce ' exp{-(y+ + y_ - y/y + y_)}t dy/coshy. 



(6.4) 
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In the second stage, different fluid elements with different rapidities switch to conic flight at different 
time t m (y). The rapidity is frozen after t > t m (y). The entropy element after freezeout needs to be 
evaluated at the switching time t = t m (y) 



dS = ce 



3/4 



exp{-(y + + y_ - y/y + y_)}- 



dy (6.5) 

t=t m (y) 



coshy 

To evaluate the square-bracketed quantity at t = t m (y), we obtain from Eq. (13.71) and (I6.2p that 

t p ±y 

* y± = (6-6) 



Therefore, we have 



which gives 



We note that 



A cosh y 



e v±\ = tjjj)_j^ = 2a e±y 

\t=t m (y) £ cos hy A V ' ^ 



l/ ± |« m(v) =ln(2fl/A)±y. (6.8) 



ln(2a/A) =y b = L + ln2, (6.9) 
where yb is the beam rapidity in the center-of-mass system, 

y b = cosh -1 (y/s NN /2m p ) = In (y/s NN /m p ) . (6.10) 
The entropy element of Eq. (16. 5p is therefore 



dS = ce 3 /i 2aexp{-2y b + yjy* - y 2 }dy (6.11) 

As the entropy is proportional to the number of particles, we obtain the rapidity distribution of the 
particle number 



^j- oc ewi^vi 2 -V 2 }- ( 6 - 12 ) 

which differs from Landau's rapidity distribution of Eq. (11.51) . with the beam rapidity y b = L + In 2 
replacing Landau's L. Therefore, we find that the rapidity distribution in the center-of-mass system 
should be modified from the original Landau distribution. 

While many steps of the formulation are the same, the main difference between our formulation 
and Landau's appears to be the additional factor of 2 in Eq. (16. 7p and (I5.2p in the new formulation. 
This factor can be traced back to the factor of 4 in the ratio of 4e/3 from (e+p) on the left hand side 
of Eq. (14. 6 p and e/3 from the pressure p on the right hand side. In Landau's formulation, this factor 
of 4 is taken to be unity for an order-of-magnitude estimate of the transverse expansion. Hence, 
there can be a modification of Landau's expression for the rapidity distribution. 
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VII. COMPARISON OF LANDAU HYDRODYNAMICS WITH EXPERIMENTAL RAPIDITY 

DISTRIBUTIONS 



We would like to compare the modified distribution with the Landau distribution and experimental 



distributions for central AuAu collisions at various energies |13l . |14| . |15l |. We can evaluate a few 
quantities to get an idea of the differences. Consider collisions at y/s NN = 200 GeV. The beam 
rapidity y& is y& = 5.36, and the logarithm of the Lorentz contraction factor is L = 4.67. The 
difference between y^ and L is substantial and leads to different shapes of the rapidity distributions 
as one observes in Fig. 6. 



exp{ V y I - y 2 } 



exp{ V L 2 - y 2 } 

exp{-y 2 /2L) 




71 
71 

K*{x.4) 
JC"(x4) 
P (x3) 
P (x3) 



FIG. 6: Comparison of experimental rapidity distribution with theoretical distribution in the form of dN/dy oc expj-^/j/^ — y 2 } 
(solid curves), Landau's distribution dN/dy oc exp{^/ L 2 — y 2 } (dashed-dot curves), and the Gaussian dN/dy oc exp{— y 2 /2L} 
(dashed curves) for produced particles with different masses. Data are from [l3j] for AuAu collisions at y?s NN = 200 GeV. 



Fig. 6 gives the theoretical and experimental rapidity distributions for ir + , tt~ , K + , K~,p, and p 
fl3| . The solid curves in Fig. 6 are the results for y/s NN = 200 GeV from the modified distribution 
Eq. (16.121) with the parameter, whereas the dashed curves are the Landau distribution of Eq. 
(jl.5p . dN/dy oc exp{ \/ L 2 — y 2 } with the L parameter. The theoretical distributions for different 
types of particles have been obtained by keeping the functional forms of the distribution and fitting 
a normalization constant to match the experimental data. We observe that Landau rapidity distri- 
butions are significantly narrower than the experimental rapidity distributions, whereas the modified 
distribution Eq. (16.121) gives theoretical results that agree better with experimental data. 

As a further comparison, we show theoretical distributions calculated with the Gaussian distri- 
bution of Eq. (11.31) as the dashed curves. We find that except for the region of large rapidities, 
the Gaussian distributions is a good representation of the modified Landau distribution. The close 
similarity between the modified distribution (16.121) and the Gaussian distribution (II. 3p explains the 
puzzle mentioned in the Introduction. The Gaussian distribution and the original Landau distribu- 
tion are different distributions. Past successes of the Gaussian distribution in explaining experimental 
rapidity data 13|, uA, ll5|] arises, not because it is an approximation of the original Landau distribu- 
tion (jl.5p . but because it is in fact close to the modified Landau distribution (16.121) that derives its 
support from a careful re-examination of Landau hydrodynamics. 

We compare theoretical distributions with the ir~ rapidity distribution for collisions at various 
energies. The solid curves in Fig. 7 are the results from the modified distribution Eq. (I6.12p with 
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FIG. 7: Comparison of experimental rapidity distribution with theoretical distribution in the form of dN/dy oc exp{^/ yf — y 2 } 
(solid curves) and Landau's distribution dN/dy oc exp{^/L 2 — y 2 } (dashed curves) for produced particles at different energies. 
Data are from the compilations in [l3l |. 



the i/b parameter, whereas the dashed curves are the Landau distribution of Eq. (11.51) with the L 
parameter. The experimental data are from the compilation of [l3|. The modified distribution of 
Eq. (I6.12p appears to give a better agreement with experimental data than the Landau distribution 
when all the data points at all energies are considered. 

While the comparison of the modified distribution gives a better agreement with experimental data 
than the original Landau distribution, it should be pointed out that the agreement can be further 
improved. Landau's condition of rapidity freezeout, Eq. (15. ip . is only a close estimate that already 
goes a long way in providing a good description. It may however admit refinement beyond the first 
estimate. We can in fact use the experimental data to guide us for a better estimate of the rapidity 
freezeout condition which can be conveniently modified to be 

x{t m ) = (a, (7.1) 

where ( is of order unity but can allow small deviations from unity. In this case, the matching time 
becomes 

t m (y) = 2^/(acoshy, (7.2) 

and the rapidity distribution becomes 

dN/dy oc ex P Wiy b +(\n()/2]2-y2}. (7.3) 

The comparison in Figs. 6 and 7 for the tt^ distributions indicates that a parameter ( slightly greater 
than unity will lead to a better agreement with experiment. 



VIII. PREDICTIONS OF PARTICLE RAPIDITY DISTRIBUTION FOR LHC ENERGIES 

We can re- write the rapidity distribution of charged particles in terms of the normalized distribution 
dF/dy 

(dN ch /dy)/(N p3It /2) = [N ch /(N paxt /2)]dF/dy. (8.1) 



18 



The normalized distribution dF/dy is 



dF 
dy 



A norm exp{ \Jy\ — y 2 } for modified Landau distribution, 

Aiorm exp{ \/ L 2 — y' 2 } for Landau distribution, (8.2) 

Aiorm exp{— y 2 /2L} for Gaussian Landau distribution, 



where A n( , rm is a normalization constant such that 



J dF/dy = 1. (8.3) 



As we remarked earlier, the approximate Landau hydrodynamical model is restricted in the region of 
\y\ < yi, for the modified distribution and L < \y\ for the original Landau distribution. The neglect of 
the tail region of the distribution will not affect significantly the distribution and its normalization. 

With the knowledge of the total charged multiplicity from Fig. 1, and the shape of the rapidity 
distribution from Eq. (18.21) . we can calculate dN^/dy / (N pait /2) as a function of rapidity. Fig. 5 
gives the predicted rapidity distributions at LHC energies. For heavy-ion collisions at y/s NN = 5.5 
TeV with full stopping (£ = 1), the maximum value of dN/dy per participant pair is about 22 at 
midrapidity. For pp collisions at y/s NN = 14 TeV with £ = 0.5, the maximum dN/ dy is approximately 
24 at y = 0. The widths of the rapidity distributions are a y ~ 3. The solid curves are for the modified 
distribution, the dashed-dot curves are for the original Landau distribution, and the dashed curves 
are for the Gaussian distribution. These results provide useful theoretical data for comparison with 
future experiments. 




FIG. 8: The predicted rapidity distributions dN c h/ dy /{N pB ,rt /2) of charged particles produced in (i) pp collisions at ^/s NN = 14 
TeV with £ = 0.5, and (ii) AA collisions at y?s NN = 5.5 TeV with full stopping (£ = 1), in Landau hydrodynamics. The solid 
curves are obtained with the modified distribution, the dashed-dot curves are obtained with the original Landau distribution, 
and the dashed curves with the Gaussian distribution. 
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IX. GENERALIZATION TO NON-CENTRAL COLLISIONS 



We have considered so far the case of central collision, for which the dimensions in all directions 
before Lorentz contraction are all equal. In non-central collisions, the transverse radius a^/2 will 
depend on the azimuthal angle <fi measured relative to the X clXIS clS depicted in Fig. 1. 

We can easily generalize the Landau model to the case of non-central collisions. The first stage 
of longitudinal expansion proceeds in the same way. The initial energy density e will need to be 
adjusted as it depends on the impact parameter. The important quantity is however the temporal 
dependence of the transverse displacement. Following the same Landau arguments as in the central 
collision case, Eq. H4.6[) for the transverse displacement can be generalized to the non-central collision 
case to be 

3 K J t 2 3a /2 v ; 

where p(4>, t) is the transverse displacement at azimuthal angle <p. As a consequence, the transverse 
displacement depend on <f> and t as 

^ Aa^u ) 2 4a^cosh 2 ?/ ^ ^ 

We can set up the Landau condition for the onset of the second stage as the condition that the 
transverse displacement p(4>, t) is equal to the transverse dimension a^, 

p(4>,t m ) = tty. (9.3) 

Thus, in the case of non-central collision, the Landau condition of (15. ip is changed to 

t m (y,4>) = (at /a) x 2acoshy, (9.4) 

which shows that the rapidity freezeout time t m is reduced by the factor (a^/a), compared to the 
freezeout time for central collisions. For a fixed azimuthal angle 0, fluid elements go from the first 
stage motion to the second stage motion at different times. The matching time t m (the rapidity 
freezeout time) occurs at different instances for different rapidities. Following the same argument as 
before, Eq. (16.81) for the non-central collision case becomes 

y±\ t=tmM) = ln((ty/a) + ln(2a/A 6 ) ± y, (9.5) 

where the longitudinal thickness of the initial slab depends on the impact parameter b. As a 
consequence, the rapidity distribution for this non-central collision is 

dN 



oc exp{ A /ln(2a/A b ) + ln(a^/a)] 2 - y 2 }. (9.6) 



X. CONCLUSIONS AND DISCUSSIONS 



In many problems in high-energy collisions such as in the description of the interaction of the 
jet or quarkonium with the produced dense matter, it is desirable to have a realistic but simple 
description of the evolution of the produced medium. Landau hydrodynamics furnishes such a tool 
for this purpose. 
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Recent successes of Landau hydrodynamics [13|, [14|, [15[ in explaining the rapidity distribution, 
total particle multiplicities, and limiting fragmentation indicate that it contains promising degrees 
of freedom. Questions are however raised concerning the use of pseudorapidity or rapidity variables, 
the approximate Gaussian form or the square-root exponential form of the rapidity distribution, and 
the values of the parameters in the rapidity distribution. 

We start with the rapidity variable from the outset so that we do not need to worry about the 
question of the rapidity or the pseudorapidity variable. We follow the formulation of the Landau 
hydrodynamics by keeping careful track of the numerical constants that enter into the derivation. 
We confirm Landau's central results except that the approximate rapidity distribution obtained 
by Landau needs to be modified, when all numerical factors are carefully tracked. In particu- 
lar, the rapidity distribution in the center-of-mass system should be more appropriately given as 
dN/dy oc exp{^y 2 — y 2 }, where y^ is the beam nucleon rapidity, instead of the Landau original 
result of dN/dy( Landau) oc exp{^/L 2 — y 2 }. The modified distribution leads to a better description 
of the experimental data, thereby supports the approximate validity of Landau hydrodynamics as 
a description of the evolution of the produced bulk matter. Phenomenological fine-tuning can be 
further introduced as in Eqs. (17. 2p and (17. 3p to relate the rapidity freezeout time t m with the rapidity 
and transverse radius. 

The modified distribution differs only slightly from the Gaussian distribution dN/ dy( Gaussian) oc 
expj— jy 2 /2L|, that has been used successfully and extensively in the literature [13 



171 ] . This explains the puzzle we mention in the Introduction. Even though the Gaussian Landau 



distribution (11.31) is conceived as an approximate representation of the original Landau distribution 
(ll.Sp for the region of small rapidity with \y\ <C L, it differs from the original Landau distribution 
in other rapidity regions. They are in fact different distributions. The Gaussian distribution has 



been successfully used to explain experimental rapidity distribution data [13|, [14|, |l5j, not because 



it is an approximation of the original Landau distribution (jl.5p . but because it is in fact a good 
representation of the modified Landau distribution (16.121) that derives its support from a careful 
re-examination of Landau hydrodynamics. Thus, there is now a firmer theoretical support for the 
Gaussian distribution (II. 3p owing to its similarity to the modified distribution of (16.121) . 

The need to modify Landau's original distribution should not come as a surprise, as the original 
distribution was intended to be qualitative. Our desire to apply it quantitatively therefore lead to a 
more stringent re-examination, with the result of the modification as we suggest. 

The quantitative successes of the modified Landau hydrodynamics make it a useful tool for many 
problems in high-energy heavy-ion collisions. We can now view with some degree of confidence 
Landau's picture for the evolution of the bulk matter. The evolution first proceeds with a one- 
dimensional longitudinal expansion with a simultaneous transverse expansion. Different parts of 
the fluid is subject to different transverse displacements. Those in the central rapidity receive the 
greatest transverse displacement and will freeze out in y at the earliest. The final distribution is 
obtained by collecting the bulk matter at different times of rapidity freezeout. 

In spite of these successes, many problems will need to be examined to make the Landau model 
an even better tool. We can here outline a few that will need our attention. The distribution so far 
deals with flow rapidity of the fluid elements, and the thermal distribution of the particles inside 
the fluid element has not been included. The folding of the thermal distribution of the particles 
will broaden the rapidity distribution and should be the subject of future investigations. Another 
improvement is to work with a curvilinear coordinate system in the transverse direction to obtain the 
transverse displacement. The will improve the description of the matching time in the transverse 
direction. One may wish to explore other forms of the freezeout condition instead of Landau's 
transverse displacement condition, to see how sensitive the results can depend on the freezeout 
condition. Finally, as the approximate solution for the one- dimensional is also available, it may also 
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be of interest to see how much improvement there can be in obtaining the matching time estimates 
which enter into the rapidity freeze-out condition. 
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